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In popular cosmological scenarios, some time beyond a redshift of 10, stars within 
protogalaxies created the first heavy elements; these systems, together perhaps with 
an early population of quasars, generated the ultraviolet radiation and mechanical 
energy that reheated and reionized the cosmos. The history of the Universe dur- 
ing and soon after these crucial formative stages is recorded in the all-pervading 
intergalactic medium (IGM), which contains most of the ordinary baryonic mate- 
rial left over from the big bang. Throughout the epoch of structure formation, the 
IGM becomes clumpy and acquires peculiar motions under the influence of gravity, 
and acts as a source for the gas that gets accreted, cools, and forms stars within 
galaxies, and as a sink for the metal enriched material, energy, and radiation which 
they eject. 
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1. Introduction 

At epochs corresponding to z ~ 1000 the intergalactic medium (IGM) is expected 
to recombine and remain neutral until sources of radiation and heat develop that 
are capable of reionizing it. The detection of transmitted flux shortward of the 
Lya wavelength in the spectra of sources at z ~ 5 implies that the hydrogen 
component of this IGM was ionized at even higher redshifts. There is some evidence 
that the double reionization of helium may have occurred later, but this is still 
controversial. It appears then that substantial sources of ultraviolet photons and 
mechanical energy were already present when the Universe was less than 7% of 
its current age, perhaps quasars and/or young star-forming galaxies: an episode of 
pre-galactic star formation may provide a possible explanation for the widespread 
existence of heavy elements (like carbon, oxygen, and silicon) in the IGM, while the 
integrated radiation emitted from quasars is likely responsible for the reionization 
of the intergalactic helium. Establishing the epoch of reionization and reheating 
is crucial for determining its impact on several key cosmological issues, from the 
role reionization plays in allowing protogalactic objects to cool and make stars, to 
determining the small-scale structure in the temperature fluctuations of the cosmic 
microwave background. Conversely, probing the reionization epoch may provide a 
means for constraining competing models for the formation of cosmic structures, 
and of detecting the onset of the first generation of stars, galaxies, and black holes 
in the Universe. 
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2. The transition from a neutral to an ionized Universe 

Popular cosmological models predict that most of the intergalactic hydrogen was 
reionized by the first generation of stars or accreting black holes at z = 7 — 15. One 
should note, however, that while numerical N-body+hydrodynamical simulations 
have convincingly shown that the IGM is expected to fragment into structures at 
early times in cold dark matter (CDM) cosmogonies (e.g. Cen et al. 1994; Zhang, 
Anninos, & Norman 1995; Hernquist et al. 1996), the same simulations are much less 
able to predict the efficiency with which the first gravitationally collapsed objects 
lit up the Universe at the end of the 'dark age' (Rees, this volume). 



(a) Photo- versus collisional ionization 

The scenario that has received the most theoretical studies is one where hydrogen 
is photoionized by the UV radiation emitted either by quasars or by stars with 
masses > 10 M©, rather than ionized by collisions with electrons heated up by, e.g. 
supernova-driven winds from early pregalactic ('Pop IIP) objects. In the former case 
a high degree of ionization requires about 13.6 x (1 + t/i lec ) eV per hydrogen atom, 
where i vec is the volume-averaged hydrogen recombination timescale, t/i Tec being 
much greater than unity already at z « 10 according to the numerical simulations 
of Gnedin & Ostriker (1997), and Gnedin (2000). Collisional ionization to a neutral 
fraction of only few parts in 10 5 requires a comparable energy input, i. c. an IGM 
temperature close to 10 5 K or about 25 eV per atom. 

Massive stars will deposit both radiative and mechanical energy into the inter- 
stellar medium of Pop III objects. A complex network of 'feedback' mechanisms is 
likely at work in these systems, as the gas in shallow potential is more easily blown 
away thereby quenching further star formation (Mac Low & Ferrara 1999), and the 
blastwaves produced by supernova explosions reheat the surrounding intergalactic 
gas and enrich it with newly formed heavy elements and dust. It is therefore dif- 
ficult to establish whether an early input of mechanical energy will actually play 
a major role in determining the thermal and ionization state of the IGM on large 
scales (Tegmark, Silk, & Evrard 1993). What can be easily shown is that, during the 
evolution of a a 'typical' stellar population, more energy is lost in ultraviolet radia- 
tion than in mechanical form. This is because in nuclear burning from zero to solar 
metallicity (Z© = 0.02), the energy radiated per baryon is 0.02 x 0.007 x tohc 2 ; 
about one third of it goes into H- ionizing photons. The same massive stars that 
dominate the UV light also explode as supernovae (SNe), returning most of the 
metals to the interstellar medium and injecting about 10 51 ergs per event in kinetic 
energy. For a Salpeter initial mass function (IMF), one has about one SN every 
150 M© of baryons that forms stars. The mass fraction in mechanical energy is 
then approximately 4 x 10^ 6 , ten times lower than the fraction released in photons 
above 1 ryd. 

The relative importance of photoionization versus shock ionization will depend, 
however, on the efficiency with which radiation and mechanical energy actually 
escape into the IGM. Consider, for example, the case of an early generation of halos 
with circular speed v c — 50 kms -1 , corresponding in top- hat spherical collapse to 
a virial temperature T v = 0.5^,m p v^./k w 10 5 ' 3 K and halo mass M — 0.1v^/GH « 
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10 9 [(1 + z)/10] _3 / 2 ft- _1 M Jf| In these systems rapid cooling by atomic hydrogen 
can take place and a significant fraction, of their total mass may be converted 

into stars over a dynamical timescale (here J7b is the baryon density parameter). 
For / = 0.05, n B h 2 = 0.02, and h = 0.5, the explosive output of 50, 000 SNe would 
inject an energy Eq ps 10 55 7 ergs. The hot gas will escape its host, shock the IGM, 
and eventually form a cosmological blast wave. If the explosion occurs at cosmic 
time t — 4 x 10 s yr, corresponding in the adopted cosmology (EdS with h — 0.5) to 
z = 9, at time At = 0.2t after the event the proper radius of the (adiabatic) shock 
is given by the standard Sedov- Taylor self-similar solution, 

Rs « ( V5 4 2/5 At 2/ 5 ^ 23 kpc (1) 

V "a / 

At this instant the shock velocity relative to the Hubble flow is 

v s « 2R s /5At ps 110 kms^ 1 , (2) 

still much higher than the escape velocity from the halo. The gas temperature 
just behind the shock front is T s — 3fim p v^/16k ps 4 x 10 5 K, more than enough 
to efficiently ionize all the incoming hydrogen. At these redshifts, it is the onset of 
Compton cooling off cosmic microwave background photons that ends the adiabatic 
stage of blast wave propagation. According to the Press-Schechter formalism, the 
comoving abundance of collapsed dark halos with mass M = 10 9 h^ 1 M© at z = 
9 is dn/dlnM ~ 5/i 3 Mpc -3 , corresponding to a mean proper distance between 
neighboring halos of ~ 40 ft -1 kpc, and to a total mass density parameter of order 
0.02. With the assumed star formation efficiency, only a small fraction, about one 
percent, of the stars seen today would have to be formed at these early epochs. 
Still, our simple analysis shows that the blast waves from such a population of 
pregalactic objects could overlap with large enough velocities to initially drive the 
intergalalactic medium to a significantly higher adiabat, T > 10 5 K, than expected 
from photoionization, and pollute the entire IGM with metal-enriched material. 
A lower density of sources - which would therefore have to originate from higher 
amplitude peaks - would suffice if the typical efficiency of star formation were larger 
than assumed here. 

Quasar-driven blast waves (Eq > 10 60 ergs) are instead quite inefficient at ioniz- 
ing the IGM, since much of the initial explosion energy is lost into the collisionless 
component (Voit 1996). They would also be too rare to fill the IGM without vio- 
lating the COBE limit on the y-distortion of the microwave background. 

(b) Cosmological H II regions 

In the following sections we will focus our attention to the photoionization of the 
IGM, i.e. we will assume that UV photons from an early generation of stars and/or 
quasars are the main source of energy for the reionization and reheating of the 
Universe, and that star formation and quasar activity occurs in collapsed galaxy 
halos. The process then begins as individual sources start to generate expanding 
H ii regions in the surrounding IGM; throughout an H n region, H is ionized and 
He is either singly or doubly ionized. As more and more sources of ultraviolet 

f This assumes an Einstein-de Sitter (EdS) Universe with Hq = 100 h kms" 1 Mpc — 1 . 
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radiation switch on, the ionized volume grows in size while the neutral phase shrinks. 
Reionization is completed when the H n regions overlap, and every point in the 
intergalactic space gets exposed for the first time to a nearly uniform Lyman- 
continuum (Lyc) background. 

When an isolated point source of ionizing radiation turns on, an ionization (I) 
front separating the H n and H i regions propagates into the neutral gas, and the 
proper volume Vi of the ionized zone grows according to the equation 

f-3^ = ^-^ (3) 

at n H tree 

(Shapiro & Giroux 1987), where N- lon is the number of ionizing photons emitted 
by the central source per unit time, nn(O) = 1.7 x 1CP 7 (fish 2 /0.02) cm~ 3 is 
today's mean hydrogen density, and all other symbols have their usual meaning. 
Most photons travel freely in the ionized gas, and are absorbed in a transition layer. 
In the case of stellar sources the I-front is quite sharp, and the degree of ionization 
changes on a short distance of the order of the mean free path for an ionizing photon. 
When t rcc <C t, the growth of the H n region is slowed down by recombinations in 
the highly inhomogencous IGM, and its evolution can be decoupled from the Hubble 
expansion. Just like in the static case, the ionized bubble then fills its time-varying 
Stromgrcn sphere after a few recombination timcscalcs, 

Vi = — = (1 - e t ' t ^). (4) 

While the volume that is ionized depends on the luminosity of the central source, 
the time it takes to produce an ionization-bounded region is only a function of t Iec . 

In the presence of a population of ionizing sources, the transition from a neutral 
IGM to one that is almost fully ionized can be statistically described by the evolu- 
tion with redshift of the volume filling factor (or porosity) Q of H n , He n , and 
He in regions. The radiation emitted by spatially clustered stellar-like and quasar- 
like sources - the number densities and luminosities of which may change rapidly 
as a function of redshift - coupled with absorption processes in a medium that be- 
comes more and more clumpy owing to the non-linear collapse of structures (Figure 
1), all determine the complex topology of neutral and ionized zones in the Universe 
(Gnedin 2000; Ciardi et al. 2000; Abel, Norman, & Madau 1999). When Q < 1 and 
the radiation sources are randomly distributed, the ionized regions are spatially 
isolated, every UV photon is absorbed somewhere in the IGM, and the UV radi- 
ation field is highly inhomogeneous. As Q grows, the crossing of ionization fronts 
becomes more and more common, until percolation occurs at Q = 1. 

Since the mean free path of Lyc radiation is always much smaller than the hori- 
zon (this is also true after 'overlapping' because of the residual H i still present 
in the Lya forest clouds and the Lyman-limit systems), the filling factor of cos- 
mological H n regions is equal at any given time t to the total number of ioniz- 
ing photons emitted per hydrogen atom by all radiation sources present at ear- 
lier epochs, J n mn dt' /Sh, minus the total number of radiative recombinations per 
atom, L Qdt' /t rcc . This statement reflects the simple fact that every ultraviolet 
photon that is emitted is either absorbed by a newly ionized hydrogen atom or by a 
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Figure 1. Simulating the reionization of the Universe: propagation of an ionization front 
in a 128 3 cosmological density field. A 'mini-quasar' emitting 5 x 10 53 ionizing photons 
s~ was turned on at the densest cell, in a virialized halo of total mass 10 11 M©. The box 
length is 2.4 comoving Mpc. The solid contours give the position of the front at 0.15, 0.25, 

0. 38, and 0.57 Myr after the quasar has switched on at z = 7. The underlying greyscale 
image indicates the initial H i density field. (From Abel et al. 1999.) 

recombining one. Differentiating one gets 

dQ _ n ion _ Q_ , 

dt n H ?roc 

(Madau, Haardt, & Rees 1999). It is this differential equation - and its equivalent 
for expanding helium zones - that statistically describes the transition from a neu- 
tral Universe to a fully ionized one independently, for a given UV photon emissivity 
per unit cosmological volume n lon , of the complex and possibly short-lived emis- 
sion histories of individual radiation sources, e.g. on whether their comoving space 
density is constant or varies with cosmic time. Initially, when the filling factor is 
<C 1, recombinations can be neglected and the ionized volume increases at a rate 
fixed solely by the ratio ji-ion/^-H- As time goes on and more and more Lyc pho- 
tons are emitted, radiative recombinations become important and slow down the 
growth of the ionized volume, until Q reaches unity, the recombination term satu- 
rates, and reionization is finally completed (except for the high density regions far 
from any source which are only gradually eaten away, Miralda-Escude, Haehnelt, 
& Rees 2000). In the limit of a fast recombining IGM (i ICC <C t), one can neglect 
the derivative on the left-hand side of equation (5) and derive 

Q < ^trcc, (6) 

1. e. the volume filling factor of ionized bubbles must be less (or equal) to the number 
of Lyc photons emitted per hydrogen atom in one recombination time. In other 
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words, because of radiative recombinations, only a fraction i Tec /t <C 1 of the photons 
emitted above 1 ryd is actually used to ionize new IGM material. The Universe is 
completely reionized when 

^-ion^rcc Z «H, (7) 

i.e. when emission rate of ultraviolet photons exceeds the mean rate of recombina- 
tions. 



(c) A clumpy Universe 

The simplest way to treat reionization in a inhomogeneous medium is in terms of 
a clumping factor that increases the effective gas recombination rate. In this case 
the volume-averaged recombination time is 

t rcc = [(l + 2 X )n p a B C}- 1 =0.06Gyrf^P) ' (^) " ^ C w \ (8) 



v 0.02 J V 10 / n p 

where as is the recombination coefficient to the excited states of hydrogen (at an 
assumed gas temperature of 10 4 K), \ the helium to hydrogen abundance ratio, 
and the factor C = {n p )/n p > 1 takes into account the degree of dumpiness of 
photoionized regions (hereafter do = C/10). If ionized gas with density n p filled 
uniformly a fraction 1/C of the available volume, the rest being empty space, the 
mean square density would be {n p ) = n p /C = fi p C. More in general, if f m is the 
fraction of baryonic mass in photoionized gas at an overdensity S relative to the 
mean, and the remaining (underdense) medium is distributed uniformly, then the 
fractional volume occupied by the denser component is 

fv = fm/S, (9) 
the density of the diffuse component is 

n P : _ f , (10) 

1 Jv 

and the recombination rate is larger than that of a homogeneous Universe by the 
factor 

c = /^+ (1 ~ / 7 )2 ( n ) 

t Jv 

(e.g. Chiu & Ostrikcr 1999; Valageas & Silk 1999). It is difficult to estimate the 
clumping factor accurately. According to hydrodynamics simulations of structure 
formation in the IGM (within the framework of CDM-dominated cosmologies), 
Lya forest clouds with moderate overdensities, 5 < S < 10, occupy a fraction of 
the available volume which is too small for them to dominate the clumping at high 
redshifts (e.g. Zhang et al. 1998; Theuns et al. 1998). In hierarchical clustering 
models, it is the virialized gas (with S w 180 if one ignores the slope of the density 
profile) in dark matter halos with temperatures < 10 4 K (masses M < 10 7 h~ Y Mq) 
which will plausibly boost the recombination rate by large factors as soon as the 
collapsed mass fraction exceeds 0.5%. Halos or halo cores which are dense and 
thick enough to be self-shielded from UV radiation will stay neutral and will not 
contribute to the recombination rate. This is also true of gas in more massive 
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halos, which will be virialized to higher temperatures and ionized by collisions with 
thermal electrons. With a large comoving space density at z = 9 of dn/dh\M ~ 
1000 h 3 Mpc~ 3 , corresponding to a mean proper distance of only ~ 6 ft. -1 kpc, and 
to a mass fraction of 0.04, halos with T v « 10 4 K will contribute significantly, 
f m S m 7, to the clumping. Recent calculations by Benson ct al. (2000), which 
instead include all halos with T v > 10 4 K and adopt an isothermal density profile 
with a flat core, give C w 30 already at z = 9. Because of finite resolution effects, 
numerical simulations may underestimate clumping: in those of Gnedin & Ostrikcr 
(1997), for example, C rises above unity at z < 20, and grows to C ~ 10 (40) at 
z « 9 (5). 

It is important to note that the use of the volume-averaged clumping factor in 
the recombination timescale is only justified when the size of the H n regions is much 
larger compared to the scale of the clumping, so that the effect of many halos within 
the ionized volume can be averaged over. This will be a good approximation either 
at late epochs, when the H n zones have had time to grow (or when overlapping 
ionized regions from an ensemble of sources are able to proper sample the small-scale 
density fluctuations), or at earlier epochs if the ionized bubbles are produced by 
more luminous sources like quasars or the stars within halos collapsing from high-<r 
peaks. As mentioned above, the mean free path between halos having T v w 10 4 K 
is A ~ 6 h~ x kpc at z = 9, but their virial radius is only r v ~ 0.4 h^ 1 kpc. It is only 
on scales greater than A 3 /r^ w 2hT x Mpc that the clumping can then be averaged 
over, and the covering factor of halos within the Stromgren sphere exceeds unity. 

3. Sources of UV photons 

(a) Quasars 

In recent years, several optical surveys (Warren, Hewett, & Osmer 1994; Schmidt, 
Schneider, & Gunn 1995; Kennefick, Djorgovski, & de Carvalho 1995) have consis- 
tently provided evidence for a turnover in the QSO counts. The space density of 
radio-loud quasars also appears to decrease strongly for z > 3 (Shaver et al. 1996), 
suggesting that the turnover is indeed real and not an effect on optically-selected 
QSOs induced by dust along the line of sight. The density of optically bright and 
flat-spectrum radio-loud quasars has a relatively flat maximum at 1.8 < z < 2.8, 
and declines gradually at higher redshifts (Figure 2). 

The QSO emission rate of hydrogen Lyc photons per unit comoving volume, Aq, 
is also shown in Figure 2. The procedure adopted to derive this quantity implies a 
large correction for incompleteness at high-z. With a fit to the quasar luminosity 
function (LF) which goes as <f>(L) cx L^ 13 , with [3 — 1.64 at the faint end (Pei 1995), 
the contribution to the emissivity converges rather slowly, as L 36 . At z = 4, for 
example, the blue magnitude at the break of the LF is M* ss —25.4, comparable 
or slightly fainter than the limit of current high-z QSO surveys. While a large 
fraction, about 90% at z = 4 and even higher at earlier epochs, of the ionizing 
emissivity shown in the figure is therefore produced by quasars that have not been 
actually observed, and are assumed to be present based on an extrapolation from 
lower redshifts, it is also fair to ask whether an excess of low-luminosity QSOs, 
relative to the best-fit LF, could actually boost the estimated Lyc emissivity at 
early epochs. The interest in models where the quasar LF significantly steepens 
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redshift redshift 

Figure 2. Left: comoving space density of bright QSOs as a function of redshift. The 
data points with error bars are taken from different optical and radio surveys (see Madau, 
Haardt, & Rees 1999 for details). Right: comoving emission rate of hydrogen Lyc photons 
(solid line) from QSOs, compared to the comoving rate of recombinations, nH(0)/t rC c 
(dashed lines) in a IGM with gas clumping factor C = 20, 30, 40). A EdS cosmology with 
n B h 2 = 0.02 and h = 0.5 has been assumed. Models based on photoionization by quasar 
sources appear to fall short at z ~ 5. The data points show the estimated contribution 
from Lyman-break galaxies at z ~ 3 and 4, assuming that the fraction of Lyc photons 
which escapes the dense H i layers into the galaxy halos and the IGM is 50%. 

with lookback time, and therefore predict many more QSOs at faint magnitudes 
than the extrapolation of Pei's (1995) fitting functions, stems from recent claims of a 
strong linear correlation between bulge and observed black hole masses (Magorrian 
et al. 1998), linked to the steep mass function of dark matter haloes predicted by 
hierarchical cosmogonies (e.g. Haehnelt, Natarajan, & Rees 1998; Haiman & Loeb 
1998). As discussed by Haiman, Madau, & Loeb (1999), the space density of low- 
luminosity quasars at high- z is constrained by the observed lack of red, unresolved 
faint objects in the Hubble Deep Field (HDF). Down to a 50% completeness limit 
of Vab = 29.6 (Jab = 28.6), no z > 4 quasar candidates have actually been found 
by Conti et al. (1999): by contrast, about 10 objects would be predicted by a QSO 
evolution model characterized by a steep LF with slope [3 — 2 and a comoving space 
density that remains constant above z = 2.5 instead of dropping (Figure 3), and 
choosen to boost the emission rate of ultraviolet photons at z ~ 5 by a factor of 5. 
A large population of faint AGNs at high-z would still be consistent with the data 
if, at these faint magnitude levels and high image resolution, the host galaxies of 
active nuclei could actually be resolved by the Hubble Space Telescope. 

(b) Star-forming galaxies 

Galaxies with ongoing star- formation are another obvious source of Lyc photons. 
The recent progress in our understanding of faint galaxy data made possible by the 
identification of star-forming galaxies at 2 < z < 4 in ground-based surveys and 
in the HDF has provided new clues to the long-standing issue of whether galaxies 
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Figure 3. Theoretical number-magnitude relation of quasars in the redshift range 
4.0 < z < 5.5. The solid line shows the prediction for a 'standard' QSO model, one 
in which the faint end of the QSO luminosity function has slope j3 = 1.64 and tracks the 
turnover observed in the space density of bright quasars at z > 3. The dashed line shows 
the increased number of sources expected in the case of a steeper, (3 — 2, luminosity func- 
tion and a comoving space density that stays constant above z = 2.5. The latter evolution 
scenario provides, within the errors, enough UV photons to keep the Universe ionized at 
z ~ 5, but appears to be inconsistent with the lack of red, faint stellar objects observed 
in the Hubble Deep Field. 

at high redshifts can provide a significant contribution to the ionizing background 
flux. Since the rest-frame UV continuum at 1500 A (redshifted into the visible 
band for a source at z ~ 3) is dominated by the same short-lived, massive stars 
which are responsible for the emission of photons shortward of the Lyman edge, 
the needed conversion factor, about one Lyc photon every 10 photons at 1500 A, is 
fairly insensitive to the assumed IMF and is independent of the galaxy history for 
t > 10 7 ' 3 yr. 

Composite ultraviolet luminosity functions of Lyman-break galaxies (LBG) at 
z k 3 and z w 4 have been recently derived by Steidel et al. (1999). They are based 
on a large catalog of spectroscopically and photometrically selected galaxies from 
the ground-based and HDF samples, and span about a factor of 40 in luminosity 
from the faint to the bright end. Integrating these LF over all luminosities L > 
0.1L*, and using the conversion L(1500)/L(912) w 6 valid for a Salpeter mass 
function and constant star formation rate, we derive for the comoving emissivities 
at 1 ryd the values of 9 ± 2 x 10 25 h ergss" 1 Hz" 1 Mpc~ 3 at z w 3, and 7 ± 2 x 
10 25 h ergss -1 Hz -1 Mpc~ 3 at z « 4, about 4 times higher than the estimated 
quasar contribution at z = 3. These numbers do not include any correction for 
local H i absorption (since the color excess -E912-1500 is expected to be small, dust 
exinction can probably be neglected in correcting from observed rest-frame far- 
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UV to the Lyman edge). The data points plotted in Figure 2 assumes a value 
of /esc = 0.5 for the unknown fraction of Lyc photons which escapes the dense 
sites of star formation (not included in our clumping factor) into the halos and 
the intergalalactic space. Note that, at z = 3, Lyman-break galaxies radiate more 
ionizing photons than QSOs for / esc > 25%. 



4. Implications 

(a) First light 

We have seen in the previous sections that, in the approximation the clumping can 
be averaged over, only the photons emitted within one recombination timescale 
can actually be used to ionize new material. As t ICC C t at high redshifts, it is 
possible to compute using equation (7) a critical value for the photon emission rate 
per unit cosmological comoving volume at a given epoch, Af c , independently of the 
(unknown) previous emission history of the Universe: only rates above this value 
will provide enough UV photons to keep the IGM ionized at that epoch. Equation 
(7) can then be rewritten as 

*m-^-<^"^>^)"(w)'- <12 > 

The uncertainty on this value is difficult to estimate, as it depends on the clump- 
ing factor and the nucleosynthesis constrained baryon density. It is interesting to 
convert this rate into a 'minimum' star formation rate per unit (comoving) volume, 
p* (for n B h 2 = 0.02): 

• AXl0 " 3 -' (0.12M oyt ->Mpc-3)(MWi±£) S . (13 , 



Jesc \ ,/csc / \ J-U 

(The conversion factor can be understood by noting that, for each 1 M & of stars 
formed, 8% goes into massive stars with M > 2OM that dominate the Lyc lumi- 
nosity of a stellar population. At the end of the C-burning phase, roughly half of 
the initial mass is converted into helium and carbon, with a mass fraction released 
as radiation of 0.007. About 25% of the energy radiated away goes into Lyc photons 
of mean energy 20 eV. For each 1 Mq of stars formed every year, we then expect 
0.08 x 0.5 x 0.007 x 0.25 x M c 2 /2OeVyr - lO^phots- 1 to be emitted shortward 
of 1 ryd.) 

Taken at face value, equations (12) and (13) have perhaps a surprising implica- 
tion. In a inhomogcneous Universe, early reionization at z ~ 9 requires an ionizing 
emissivity which is comparable or larger than that radiated by QSOs at the peak 
of their activity, z w 3. In a similar manner, photoionization by massive stars can 
only play a role if the star formation density at this epoch were significantly larger 
than the value directly 'observed' (i.e. uncorrected for dust reddening) at z = 2 
(Madau, Pozzetti, & Dickinson 1998). 



(b) Delayed He n reionization 

Because of its higher ionization potential and the steep spectra of UV radia- 
tion sources, the most abundant (by a factor ~ 100) absorbing ion in the post- 
rcionization Universe is not H i but He n . The importance of intergalactic helium 
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Figure 4. The evolution of the He m filling factor as a function of redshift in a inhomoge- 
neous Universe where photoionization is dominated by QSOs turning over at z > 3. From 
right to left, the three curves assume a constant clumping factor of C = 10, 20, and 30. 
The QSO photon spectrum is assumed to vary as v~ 2 S shortward of the hydrogen Lyman 
edge. Note how the ionization of He n is never completed before z — 3 in models with 
C > 10. 

in the context of this study stems from the possibility of detecting the effect of 
'incomplete' He n reionization in the spectra of z ~ 3 quasars as, depending on 
the dumpiness of the IGM (Madau & Meiksin 1994), the photoionization of singly 
ionized helium may be delayed until much later than for H i . 

Since H i and He i do not absorb a significant fraction of hv > 54.4 eV photons, 
the problem of He n reionization can be decoupled from that of other ionizations, 
and the equivalent of equation (5) for expanding He in regions becomes 

dQ _ n ion4 Q 

OX riHo IHoIII 

where n; otl 4 now includes only photons above 4 ryd, and tnoin is 6.5 times shorter 
than the hydrogen recombination timescale if ionized hydrogen and doubly ionized 
helium have similar clumping factors .| It is interesting to note that, if the intrinsic 
photon spectrum of ionizing sources has slope h(y) oc z^~ 2 ' 8 , the first terms on 

t This last assumption appears, however, rather dubious: the reason is that self-shielding of 
He II Lyc radiation occurs at much lower hydrogen columns than self-shielding of photons at 1 ryd 
(by about a factor of S/2, where the spectral 'softness' S is the the ratio of the radiation flux at 
the hydrogen Lyman edge to the flux at 4 ryd), and self-shielded gas will remain neutral and not 
add to the recombination rate. Ionized hydrogen may then be more clumpy than doubly ionized 
helium. 
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the right-hand side of equations (5) and (14) are actually equal, and a significant 
delay between the complete overlapping of H n and He in regions can only arise 
if recombinations are important. This effect is illustrated in Figure 4, where the 
expected evolution of the He in filling factor (obtained by numerical integration 
of eq. 14) is plotted for a QSO-photoionization model with a source decline at 
high rcdshifts: He n reionization is never completed before z — 3 in models with 
C > 10. A significant contribution to the UV background at 4 ryd from massive 
stars, which could push the helium reionization epoch to higher redshifts, has been 
traditionally ruled out on the basis that the ratio between the number of He n and 
H i Lyc photons emitted from low-metallicity starbursts is only about two percent 
(Leitherer & Hcckman 1995), five times smaller than in typical QSO spectra. It has 
been recently pointed out by Tumlinson & Shull (2000), however, that metal-free 
stars exhibit higher effective temperatures and dramatically harder stellar spectra, 
particularly in the He n continuum. This enhanced He-ionizing capabilities of Pop 
III stars could have interesting implications for reionization. 

To date, various studies of the Hell Lya forest in the spectra of distant QSOs 
(Hogan, Anderson, & Rugcrs 1997; Rcimers et al. 1997; Heap et al. 2000) have 
revealed patchy absorption with low He n opacity 'voids' alternating several Mpc 
sized regions with vanishing flux. These observations suggest that helium absorption 
does not increase smoothly with lookback time, but rather in the abrupt manner 
expected in the final stages of inhomogencous reionization by quasar sources. Ra- 
diative transfer effects during He n reionization could affect the thermal history 
of the IGM (Abel & Haehnelt 1999, Efstathiou, this volume). Here it is impor- 
tant to remark that, while delayed He n reionization in a clumpy Universe appears 
to be naturally linked to the observed decline in the space density of quasars be- 
yond z ~ 3, the complete overlapping of He in regions occurs instead much earlier 
(z ^> 5) in models that predict many more faint QSOs at high redshifts (Haiman 
& Loeb 1998). 

I would like to thank my collaborators, T. Abel, F. Haardt, Z. Haiman, and M. 
Rees, for many useful discussions on the topics discussed here. 
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